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Graphene exhibits zero-gap massless Dirac fermions and zero density of states at = 0.^ 
These particles form locahzed states called edge states on a finite-width strip with zigzag 
edges at E = 0.^ Naively thinking, one may expect that a current is concentrated at the edge, 
but Muhoz-Rojas et al. noticed that it is homogeneously distributed,'^ and Zarbo and Nikolic^ 
numerically obtained a result that the current density is maximum at the center of the strip. 
We rigorously derive an expression for the current density and clarify the reason why the 
current density of the edge state has a maximum at the center. 

The expectation value of local current density for a state is defined as 

2et 

{Mj(x,i),(x',i')\'4^) = -^I™ [K'{x')ai{x)] , (1) 

where j(x,i),{x',i') is the current density operator of the bond from site {x,i) to site {x',i'), 
ai{x) and bii{x') are the corresponding amplitudes of the wave function at these sites as 
shown in Fig. 1, e is the charge, and t is the nearest-neighbor hopping parameter. Because 
of the translational symmetry along the x-axis, the amplitudes of the wave function at the 
A-sublattice and B-sublattice with wave number k are respectively defined as 

ai{x) = Aiex.p{ikx), (2) 

bi{x) = Biexp{ikx). (3) 

The Schrodinger equations in the tight-binding model with nearest-neighbor hopping are 
written as 

KA, + tA+i = EBi, (4) 

tB, + KB,+i = EAi+i. (5) 

Here, K = 2tcos {ka/2), E is the energy, and a is the lattice constant.^ In the region 27r/3a < 
< vr/a, a solution with ~ exists, which is the edge state. ^ In this region, the wave 
amplitudes at the A- and B-sublattice are approximately written as 

Ai = All -2cos — j , il<i<w) (6) 

/ ka\ 

Bi = ±Ai ( -2 cos Y j . {l<i<w) (7) 



1/5 



J. Phys. Soc. Jpn. 



Short Note 



In these equations w defines the width of the strip. The wave function becomes exact in the 
hmit of either w ^ oo or k ^ vr/a. 




bi+i {x - f ) bi+i {x + f ) 



Fig. 1. (Color online) Definition of wave function. Part of the strip is sliown. fn this structure am- 
plitude of the upper (lower) edge state is mainly on the A-sublattice (B-sublattice). 



In this paper, we derive a relation for the local current density in the x-direction. The 

local current density flowing through the bond between Oj (x — a/2) and bi{x) is written as 

2et . „ ka 
3i = ^AiBism—. (8) 

Using eqs.(4), (5), and (8), we obtain the relation between ji and jj+i as 

ji+i - ji = - Bf) tan ^. (9) 

At i = 1, namely, at the upper zigzag edge, 

ji = — A^tan— . (10) 

The y-dependence of the local current density is calculated from the initial condition [eq.(lO)] 
and the recurrence equation [eq.(9)]. In the edge state, where \E{k)\ < E{2iT/3a) = 
2tsm[TT/2{2w + 1)],6 the wave function is localized around the edges and is approximately 
given by eqs.(6) and (7).'^ Notice that the wave function at the A-sublattice is localized at the 
upper edge and that that at the B-sublattice is localized at the lower edge as shown in Fig. 
2. Therefore, A'^^-^ — Bf > in the upper half of the strip, and the local current density ji 
increases as i increases. At the center of the zigzag strip, ji is maximum since Af+i — B^^ = 0. 

When If ^ 1 or A; ~ 7r/a, the coupling of edges becomes small, the wave functions given 
by eqs.(6) and (7) become almost exact, and E approaches zero. In such a situation eq.(9) 
implies that ~ jj, namely, the current flows almost uniformly. This is also implied by 
eqs.(6) and (7). By substituting eqs.(6) and (7) into eq.(8), ji becomes independent of i, and 
is written as 

2et ( ka\^^^ ,9 . ka 



jE^o = ( -2 cos — ) Ai sin — . (11) 
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Fig. 2. (Color online) Wave function of the edge state. The size of each circle shows the amplitude of 
the wave function at that site, the open red (closed blue) circles reperesent the amplitude of the 
wave function on the A-sublattice (B-sublattice). 



The energy of the edge state in such a situation is obtained from eqs.(4) and (5) to be® 

E{k) ~ ±t 



ka 

-2 cos — 



1 



ka\' 
2COS-J 



(12) 



The expression for the current is consistent with that obtained from the group velocity 



v{k) 



dE{k) at 
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(13) 



at this energy.^ When \E\ — > 0, the magnitude of the current for a state becomes small, 
as implied by eq.(ll). However, the conductance is still 2e^//i, because the smallness of the 
current is compensated for by the large density of states. 

To verify our analytical discussion, we numerically calculated the current density from 
eq.(l). Using Zarbo and Nikolic's parameter {E = O.Olt, if = 10), we reproduced their resulf^ 
as shown in Fig. 3(a). With decreasing energy, the current flows almost uniformly in the 
y-direction as shown in Fig. 3(b). These results are consistent with those derived from eq.(9). 

In summary, a finite-width strip with zigzag edges has edge states with energy |i?(A;)| < 
E{27r/3a), and electrons are strongly localized around the edges. However, the current flows 
almost uniformly at the center and decreases towards the edges. This behavior is understood 
by the expression for the local current density in the x-direction. This uniform current is 
realized as a result of three factors. Firstly, the wave function of the edge states decays 
exponentially towards the other edge. Secondly, graphene consists of an A-sublattice and B- 
sublattice. Thirdly, the local current density is calculated by the multiplication of amplitudes 
on both sublattices.^ 

In this short note we considered the current distribution of an edge state specified by wave 
number k. The interaction between electrons is not taken into account. We remark that the 
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Fig. 3. (a) Local current density ji/jmax for E = O.Olt and w = 10, where jmax is the vahie of at 
the center. This result is consistent with that of Zarbo and Nikolic.'' (b) Local current density for 
E = IQ-H and w = 8. 



possibility of magnetic ordering of the edge state has been suggested by several authors.^' 
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